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INTEGRABLE MAPS FROM GALOIS DIFFERENTIAL 

ALGEBRAS, BOREL TRANSFORMS AND 

NUMBER SEQUENCES 

PIERGIULIO TEMPESTA 



Abstract. A new class of integrable maps, obtained as lattice versions of 
polynomial dynamical systems is introduced. These systems are obtained by 
means of a discretization procedure that preserves several analytic and alge- 
braic properties of a given differential equation, in particular symmetries and 
integrability 1401 . Our approach is based on the properties of a suitable Ga- 
lois differential algebra, that we shall call a Rota algebra. A formulation of 
the procedure in terms of category theory is proposed. In order to make the 
lattice dynamics confined, a Borel regularization is also adopted. As a byprod- 
uct of the theory, a connection between number sequences and integrability is 
discussed. 
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1. Introduction 



The analysis of evolution equations on a discrete background is a very active 
research area, due to their ubiquitousness both in mathematics and in theoretical 
physics. The formulation of quantum physics on a lattice is motivated, for in- 
stance, by the need for regularizing divergencies in field theory [l^, [2T], [28] and 
in several scenarios of quantum gravity [3], [17], [42], [36], [13]. This would imply 
the discreteness of space-time geometry. Discrete versions of nonrelativistic quan- 
tum mechanics have also been proposed (see, for instance, |12j . |45j and references 
therein). 

Integrable discrete systems have been widely investigated as well. For several 
respects, they seem to be more fundamental objects than the continuous ones. For 
this reason, many efforts have been devoted to the construction of discrete systems 
possessing an algebro-geometric structure reminiscent of that of continuous models. 
In particular, a challenging issue is to discretize nonlinear ordinary and partial 
differential equations in such a way that symmetry and integrability properties 
be preserved (see the recent interesting monograph |38j for the treatment of the 
Hamiltonian point of view). Among the classical examples of discrete integrable 
systems are the Toda systems and the Ablowitz-Ladik hierarchies [I]-[2]. 

A rich literature exists on discrete differential geometry and related algebraic 
aspects [4], [5], [30], [18], [11], [M], [19], [20], [26], [29], [ST], [32], [37]. Frobenius 
manifolds are also relevant in the discussion of generalized Toda systems [8], [9]. 
The problem of physically consistent discretizations has been considered in the 
context of field theories and Hamiltonian gravity, for instance in J13l , [15] , |36j . 

In the paper [ID], a new approach, offering a possible solution to this problem 
has been proposed. Indeed, by means of this approach one can preserve both Lie 
symmetries and Lax pairs associated to a given nonlinear PDF. Also, it allows to 
constructs new hierarchies of integrable discrete equations, representing discrete 
versions of the Gelfand-Dikii hierarchies. In this context, the notion of integra- 
bility is equivalent to the requirement that a large class of exact solutions can be 
analytically constructed, starting from those of the continuous models we wish to 
discretize. In particular, analytic solutions of the continuous systems are converted 
into exact solutions of the corresponding systems on the lattice. 

In the present work, by extending the results of j40^ , we focus on the discretiza- 
tion of continuous dynamical systems. Precisely, the problem we address is that of 
the discretization of polynomial vector fields of the form: 



(1) —z = aiqz^ + aN_iz^ -^ + . . . + fliz + ao. 



with iV e N, z = z{t) e M, oo, . . . , Oat e R. 

To this aim, a suitable Galois differential algebra is introduced. It is a pair 
{J-Ci Q), where {Fc, +,-,*) is an associative algebra of formal power series over a 
lattice £, endowed with an associative and commutative product "*", and Q is a 
delta operator that acts as a derivation with respect to this product, by preserving 
the Leibnitz rule (see also [13], [B]). A sequence of basic polynomials is uniquely 
associated with Q. An analogous product also has appeared, in a different context, 
in the theory of linear operators acting on spaces of polynomials [TB] . 
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The theoretical scheme we propose has a natural formulation in the language of 
Category Theory. We introduce the category Qal{F, Q) of Galois differential alge- 
bras and that of abstract polynomial dynamical systems £. The correspondence 
among continuous and discrete equations is expressed in terms of morphisms of 
algebras. A functor F : Qal -^ £ la defined. Essentially, it enables to map differ- 
ential equations into difference equations by preserving the underlying differential 
structure: the symmetry and integrability properties of a continuous model are 
naturally inherited by the discrete ones associated and solutions are mapped into 
solutions. This construction is made explicit in the case when the delta opera- 
tor Q is a forward difference operator (Theorem [Til) . Consequently, the algebraic 
structure underlying the dynamical models involved becomes transparent. 

Other approaches existing in the literature deal with the discretization of dif- 
ferential equations on suitable nonlinear lattices, adapted to the symmetries of the 
problem (see, for instance, the review |26|). Instead, in this work we introduce a 
new family of integrable maps, defined on a regular lattice of equally spaced points. 
This kind of lattice is usually more convenient for physical applications, and is 
especially suitable for the implementation of numerical integrating algorithms. 

The theory, being defined on a suitable associative algebra !F, is intrinsically 
a nonlocal one: the value of a dynamical observable on the lattice depends on 
different points of the lattice. This is actually a very common feature observed in 
the theory of discrete integrable models. When the dependence is on an infinite 
number of points, the approach still holds, but in a formal sense. A regularization 
procedure should be adopted, in order to give a physical content to this case. This 
resembles very much what happens in nonlocal theories when dealing with operator 
product expansions: all infinite contributions of every operator should be taken into 
account. 

However, in several cases, the value of an observable will depend on a finite 
number of lattice points, and the discretization will be said to be effective. With 
this terminology, we mean that the solutions obtained possess the general form 

n 

(2) 2„ = ^afe/(fc), 

fe=0 

where a^ G M or C and f{k) is a function of the point k on the lattice. 

An interesting byproduct of our approach is a connection between integrability 
and the theory of number sequences. Suppose we have a recurrence relation involv- 
ing several points on a lattice. This recurrence relation defines a discrete dynamical 
system in an auxiliary space of variables, that we assume to be integrable. Once we 
obtain a solution (a particular one) of the dynamical system, usually involving a 
number sequence, we have constructed a solution of the original recurrence relation. 
A different connection between number theory and difference operators, related to 
the theory of formal groups, has been established in [39], [41]. 

The necessary requirement to make it possible an effective discretization proce- 
dure is that the point of the lattice C we work over should correspond to the set 
of zeroes of the basic polynomials associated with Q. Therefore, for higher order 
delta operators, nonlinear lattices should be used. A procedure to compute the 
basic polynomials has been proposed in [10] , [21] • In full generality, the problem of 
the determination of their zeroes in a closed form is essentially open [23] . 
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The theory developed here can be easily adapted to the study of physically rele- 
vant models, especially nonlinear lattice field theories, as for instance the Liouville 
theory. It can also be developed to construct in a novel way quantum mechani- 
cal models on the lattice. The problem of a symmetry preserving discretization 
of the Schrodinger equation has been considered, for instance, in [24], [25]. The 
use of Galois differential algebras would provide a novel manner to solve the prob- 
lem, possibly more suitable for the applications, since the correspondence with the 
continuum case would be quite direct. 

Another interesting open problem is to extend the proposed technique to the dis- 
cretization of isochronous dynamical systems, in the spirit of the general framework 
proposed in [7]. 

The structure of the paper is the following. In Section II, an introduction to 
the algebraic techniques relevant for our discretization procedure is proposed. In 
Section III, the notions of Rota algebra and the categories Qal{T, Q) and £ are 
introduced and the main Theorem 1111 proved . In Section IV, some explicit examples 
of dynamical systems obtained according to the previous theory are constructed and 
some number theoretical aspects related to our approach are presented. 

2. Integrability-preserving discretizations of differential 
equations: a general framework 

2.1. Basic definitions. In this Section we review some basic notions concerning 
the algebraic theory of polynomial sequences and of finite difference operators, much 
in the spirit of the monographs |33]-[35j. 

Let us denote by {pn {x)}nefi-, n = 0, 1, 2, . . . , a sequence of polynomials in x of 
order n. Let J^ denote the algebra of formal power series in one variable x G K, 
where K is a field of characteristic zero, endowed with the operations of sum and 
multiplication of series, and product of a series by a scalar. An element of J-' is 
expressed by a formal power series of the form 

oo 

(3) /(x)^^^^ 

fc=0 

Let T be the shift operator, whose action is given by Tf (x) = f {x + a). The 
operator T can also be represented in terms of a differential operator as T = e"'^ , 
where D denotes the standard derivative. 

Definition 1. A linear operator S is said to be shift-invariant if commutes with 
T. A shift-invariant operator Q is called a delta operator if Qx — const ^ 0. 

We shall denote by V the set of delta operators. 

Definition 2. A polynomial sequence Pn (x) is called a sequence of basic polyno- 
mials for a delta operator Q if it satisfies the following conditions: 

1) po{x) = 1; 

2) Pn (0) = for all n > 0; 

3) Qpn (x) = npn-i {x) . 

Every delta operator Q has a unique sequence of associated basic polynomials 
[35] . Besides the standard derivative operator dx, a general class of difference 
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operators can be defined by [24] 

^ m 

(4) Ap = -y^akT'', l,me'Z, I < m, m-'l=p, 

k=l 

where a is the constant lattice spacing and a^ are constants such that 

m m 

(5) ^afc=0, ^fcafc=c. 

k=l k=l 

and a„j 7^ 0, a; 7^ 0. We choose c = 1, to reproduce the derivative dj. in the 
continuum limit. A difference operator of the form Q, which satisfies the equations 
([51), is said to be a delta operator of order p, if it approximates the continuous 
derivative up to terms of order a^ . 

As eq. (U]) involves m — I + 1 constants at, subject to just the two conditions 
([5|), we can fix all constants ak by choosing m — I — \ further conditions. 

2.2. Associative algebras and nonlocal functional products. Let Q be a 

delta operator, and {pn{x)}nen be its associated basic sequence. Let us denote by 
V the space of polynomials in one variable x. Since the basic polynomials Pn{x) 
for every Q provide a basis of J^ , any f & T can be expanded into a formal series 
of the form 

00 

(6) /(x) = ^a„p„(a;). 

Let £ be a lattice of equally spaced points on the line, indexed by an integer 
variable x. We shall denote by J^£ the vector space of the formal power series 
defined on C. We can endow the space F (and consequently Fc) with the structure 
of an algebra, by introducing a suitable product related with the delta operator Q 
involved in the discretization considered. Precisely, in [40], the following result has 
been proved. 

Theorem 3. For any delta operator Q, whose associated sequence is {pn{x)}nefi, 
the product 

* -.V xV — >M. 
defined as 

Pn{x) *Pmix) ■.= Pn+niix) 

and extended by linearity on the lattice, is associative, commutative and satisfies 
the Leibnitz rule: 

(7) Q {fix) * g{x)) = Q{f{x)) * g{x) + f{x) * Q{g{x)), 
where f{x),g{x) £ F. 

In the following, with a slight abuse of notation we will use the symbol "*g" 
whenever we wish to emphasize the dependence of the "*" product on the choice 
ofQ. 

One of the main features of the * product is the fact that it is nonlocal. Indeed, 
the product f * g does depend on the values of / and g at distinct points of the 
lattice, with the exception of the standard point-wise product of functions, that 
corresponds to the choice Q = dx- The simplest version of this product, i.e. that 
associated with the forward difference operator A, has been proposed in [43] . 
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In general, the spaces T are infinite-dimensional. Nevertheless, the expansion 
([6]) can also truncate in some cases, that are of special interest for the applications. 

3. Category theory, differential equations and their discretization 

The discretization approach we wish to propose can be formulated in a natural 
way in the context of Category Theory. To this aim, we need to propose several 
new definitions. 

3.1. Rota algebras. We introduce first the notion of a Rota algebra, as the natural 
Galois differential algebra over which the discretization procedure we are interested 
in is carried out. To motivate this definition, crucial for the subsequent discussion, 
we shall recall first that of Rota correspondence, that has been discussed extensively 
in the literature [35], [2l]-[25], [10] [22] ■ We will formulate it in the following way. 

Definition 4. We shall call the correspondence expressed by the following diagram 



d ^ Q 

i -I i -I 

(8) {X"}n A {pn{x)U 

the Rota (or umhral) correspondence. Here i is the isomorphism between a delta 
operator Q and its associated basic sequence {p„(x)}„gN, ip : V ^ V and (j) : V -^ 
V are applications which transform delta operators into delta operators and basic 
sequences into basic sequences respectively. 

By means of the Rota correspondence we can associate with a linear differential 
equation (defined in the usual algebra of C°° or real analytic functions, endowed 
with the point-wise multiplication of functions) a linear difference equation in such 
a way that several algebraic properties are preserved. The application of this corre- 
spondence to the Schrodinger equation, for instance, allows to discretize it in such 
a way that integrals of motion are conserved [25j . 

Nevertheless, in order to treat the more general case of nonlinear equations, one 
needs also to define their discrete versions on Rota algebras. Consequently, we shall 
formulate the Rota correspondence in the framework of category theory. 

As an immediate consequence of the Theorem [S] we have the following result. 

Proposition 5. For any choice of Q, the space {J- , +, •, *g) endowed with the sum 
of series, the multiplication by a scalar and the * product, defined as in Theorem\^ 
is an associative algebra. 

Therefore we can propose the following definition. 

Definition 6. A Rota algebra is a Galois differential algebra {J- , Q), where [J- , +, •, *g) 
is an associative algebra of formal power series, Q is a delta operator and * is the 
composition law, defined as in Theorem\^ such that Q acts as a derivation on T : 

(9) i) Qia + b)^Q{a) + Q{b), 

(10) a) Q{a*b)^ Q{a)*b + a*Q{b). 

Proposition 7. Given a delta operator Q, there exists a unique Rota algebra {J- , Q) 
associated with it. 
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Proof. It follows from the uniqueness both of the "*q" product associated with 
Q and of the sequence of polynomials {pn{x)}neN, that ensures that, once Q is 
assigned, the algebra {T, +, ■, *g) is uniquely determined. D 

3.2. The Galois category. We introduce a subcategory of the well known cate- 
gory of associative algebras 27 ; the Galois category. 

Definition 8. The Galois category, denoted by Qal{J- , Q) is the collection of all 
Rota algebras (J^, Q), with morphisms defined by 

l^Q,Q' ■■ Qal{F, Q) -^ gal{T, Q) 

(11) {^,+,;*Q)^{^',+r,*Q') 
which are closed under composition. 

The action of the morphism fJ-Q,Q' on formal power series is defined by 

(12) ^ anPnix) — > ^ amqm{x), 

n m 

where {pn{x)}neN and {qm{x)}meTi are the basic sequences associated with Q and 
Q' respectively. The property of closure under composition is trivial. 

We also introduce the new category of the abstract polynomial dynamical systems. 

Definition 9. The category of abstract polynomial dynamical systems, denoted by 
£ , is the collection of all equations of the form 

(13) Gq{Q, z, *) :— Qz — unz*^ — on-iz*^"^ — . . . + aiz — oq = 0, 

where z*^ := z * . . .*z, Oq, . . . , a^ are indeterminates over M. and A^ G N. The set 

N — times 

of correspondences 

(14) UQ^Q>:£^£, 

(15) eq(Q,z,*)^e<z(Q',z,*'), 
defines the class of morphisms of the category. 

The closure of the morphisms uq,q' under composition is easily verified. As a conse- 
quence of the previous construction, we can define a functor relating the categories 
of Galois differential algebras and abstract dynamical systems. 

Theorem 10. The application 

F : gal{F, Q) -^ E, 

(16) (J--f,-,*Q)^ 6(7(2,2,*), 

(17) MQ,S' >UQ^Q,, 

is a functor between categories. 

Proof. A direct verification shows that F preserves the composition of morphisms: 

F{t^Q",Q' O HQ',q) = F{nQn^Q,) O F{fJ,Q,^Q). 

If we denote by idg :— fiQ.Q the identity morphism, we also have 

F{idQiA))=idQiF{A)), 
wheie A (Egal{J',Q). D 
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This functor encodes all the main features of the discretization procedure we 
propose. It generalizes considerably the Rota correspondence ([8]). 

3.3. Main theorem. The following Theorem [Til represents the main result of this 
Section. As a particular instance of the action of the functor F, it guarantees the 
discretization of a differential equation into a difference equation belonging to the 
same category £, with the property that their general solutions share the same 
algebraic structure. 

Theorem 11. Given a dynamical system of the form 

(18) -—z = ajqz^ + on^iz^^^ + . . . + aiz + ao, 
at 

where z : M — > M is a C°° function, oq, . . . , a^ G K, let 

(19) Qz = aNZ*^ + aN-iz*^^-'^^ + . . . + aiz + aq 

he the abstract operator equation associated with it, defined in the Rota algebra 
{^, Q) where z G [T , +, •, *g). Let 

oo 

(20) z = Y, ht'' 

k=Q 

be a solution of (jlSp in the ring of formal power series m t G M. Then the equation 



Y.' 



f-^\ki+...kn+n 



Zn+1 — Zn — ClN } ^ f-j ; j Zk Zj. ' . . z^ ; — 

ki,...,kN=0 

+ OJV-I 2^ JT-j J^ j Zfe^Zfej • • • 2fe„-l 



fci!...fcAr_i! (n- fci - fc2 - . . . - fcAr_i) 

" , " / (_l)fel+fe2+n n\ 

+ ■•■ + «2^ E fc,!fc,. "^-"''- (n-fc,-fc2)! +"^"" + "°' 

(21) 

t/iat represents eq. (|19p /or Q = A+ on a regular lattice of points, indexed by the 
discrete variable n, admits as a solution the series 



(22) z„ = I] ^* 



A:=0 



{n-ky: 



Here we denote by y a sum ranging over all values of the indices ki such that 

J2i ^i < "-• 

Proof. The proof entails several steps. First, we prove that the equation P^ con- 
verts into the discrete map (PT|) for Q = A+. To this aim, we assume that z G C°° 
is defined on an equally spaced lattice of points £, and that t = na, with n G N 
and a G K"*" U {0}. The continuum limit is restored when cr — )■ 0. Unless specifically 
stated, in what follows we put a = 1. Also, we introduce an auxiliary space of 
variables, by means of a finite transformation on the lattice. 
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Definition 12. Letpn{t) be a basic sequence for a given delta operator Q. We call 
the transformation 

oo 

(23) z{t)^J2z„pr,{t) 

n=0 

a discrete interpolating transformation with coefficients z„ G M. 

In the specific case wlien Pn(t) are tlie lower factorial polynomials, the discrete 
transform ([23]) is finite. Indeed, we have 



(24) Pn{k) - ^ fc. 



if k < n 

It is straightforward to prove that 

nl ^ 



(25) ^n = E 



{n-iy. 



Zl 



and, for the inverse interpolating transform,, 

1 

:Zi. 



1 
(26) £„ = ^(-l^"-' 



ii{n~iy. 



1=0 

v*2 



To derive the r.h.s. of eq. pip , we start computing explicitly the product z*' = z*z. 
We get 

oo 

(z*z)„= ^ zi^zi^Pi^+i^{n) ^ 



ll,l2=0 

n n l\ I2 



^o^ok^okio ^m,-k,)lk,lil,-k,)l 



where we have introduced the kernel 

n 

(27) K„,fc,,,, =^ (-l)'i+'^^_^ . ^__^^___^. 

This expression, after some algebraic manipulations, reduces to 

(-l)'n! 2^-^^-^'' 



(28) ^nMM= E 



Putting s = I — ki — k2, and summing over s, we arrive at the final expression for 
the kernel (l27l): 



(— l)"nl 
(29) K„,fc,,fc, ^ \' — -— , n>ki + k2. 

(n- ki - k2)\ 

By generalizing the previous reasoning, we prove the first statement of Theorem 
[TTl i.e. that formula ([2T|) is the discrete analog of eq. (ITSl) on the regular lattice 
t — n. 

The second part of the Theorem, concerning the existence of the solution (|22|) 
for eq. (|2T|) . can be proved as follows. By means of the previous construction, we 
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have shown that eq. ([21]) is nothing but the image of ([TSl) under the action of the 
niorphism ug^^+ defined in E. The morphisni iiq,a+ provides the correspondence 

(30) Y.^kt" ^Y.^kPn{k). 

k k 

By means of the action of the functor F, any formal series solution z of eq. (fT9|) 
will correspond to a solution z[t) of eq. p^ and to a solution z„ of eq. (PT|) . In 
addition, on the lattice £, the sum J^k'^kPnik) truncates and converts into the 
finite sum (|22|) . 

D 

Comment 1 . The scheme offered by Theorem [10] can be used to generalize Theo- 
rem [11] to infinitely many other integrable maps, each of them defined in a specific 
Rota algebra (J^, Q). 

4. Applications: recurrences, number sequences and 

BOREL REGULARIZATION 

In this Section, a connection between integrability and the theory of number 
sequences will be established. It reposes on the fact that a recurrence relation, 
defined on a suitable space of variables, can be associated with each of the discrete 
dynamical systems previously introduced. The solutions of these dynamical models 
are strictly related to those of the recurrences, that are expressed in terms of specific 
number sequences. For completeness, in the following we will study in detail some 
particular cases of this construction. 

4.1. A new family of quadratic dynamical systems. Let z = z{t) : M ^ R a 
C°° function. Consider the dynamical system 

,31) i. = .^. 

By applying Theorem [TTJ we discretize this model in such a way that the solu- 
tions of its discrete versions conserve the same structure of the solutions of ([5T]) . 
Consequently, we associate with the system ([?T]) an abstract operator equation 

(32) Qz^z*^. 

Here z e (-7^, +, -j *q)- In order to get a discrete map, we represent the equation 
([32]) on a specific function algebra. As customarily, we denote by Zn the value of 
the function z at the point n. Eq. ([32]) becomes on the lattice 

" , ri\ki+k2+n „| 

(33) z„+i-z„ = 4 = 2^ ^-^y-j Zk.Zk. 



kilk2l ^' ^'{n-ki~k2y.' 

ki,k2=0 \ I ^J 



which is a particular instance of eq. ([21]) . 

This system possesses an interesting alter ego in the auxiliary space A of the 
variables z, defining a recurrence relation. Indeed, observe that 

oo n . 

(34) Az„ = J2 ^ilpi-i{n) - I] . "" n, (l + 1)^^+1 

fe=o ;=o ^ '' 
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n! 



and 

oo r 

(35) 4^= Yl ^hZhPh+h ^ Yl (n-lV ^ ^hzi-h- 

We conclude that the following recurrence holds 



(36) {I + l)zi+i ^ ^ zi'Zi-v 

l'=Q 

This recurrence can be solved by using the ansatz 

(37) % = zl+\ 
We deduce the solution 

n 

n 



(38) z„ = ^ 






(n-0! ° ' 
1=0 ^ ' 

which can be also expressed in the compact form as 

(39) z„ = zo •2-Fb(l,";2o)- 

Therefore, the discretization procedure we adopted is effective. Also, it preserves 
integrability. Indeed, if we write the general solution of ([5T|) in the functional space 
[T, +, •, *a), we get the corresponding solution of the discrete model ([55]). 
To prove this, observe that the general solution of eq. pip is provided by 

(40) z{t) = -^. 

1 - z^t 

We get 

(41) z{t) = zo + zlt + zlt^ + ... 

and by virtue of the correspondence <" — > Pn{t) defined by the morphism fJ.g^A+y 
taking also into account ([M]) . we obtain eq. ([55)) . 

4.2. Borel regularization. The solution ([55]) of the discrete dynamical system 
(1551) diverges as n -> oo. Instead, the solution (HO)) of the continuous model ([51]) 
has a stable point at z(oo) = 0. In order to construct a system possessing localized 
dynamics, we introduce a finite Borel-type regularization procedure. 

Definition 13. Given a series S E J-, with S ~ '^'kLo^kz'' , we call finite Borel- 
type regularization of S the following transformation of the sequence of partial sums 
{Sn}neN of S, i.e. 



n\ 

k=0 



(42) B(S„)=V&fc- 
Let 

(43) Wn = B{Zn) 



J+1 



^(n-0! " 



12 
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Consequently, the Borel regularized dynamical system associated with ([33]) is 

11, 

1 



(44) 



(71 + 1)W„+1 - W„ == ^ (-l)''^ 



+k2+n 



WkiWk2 



fci,fc2=0 



{n — ki — k2)\ 



This new dynamical system is still in correspondence with the original model pip . 
Indeed, the Borel transform of the series obtained by applying to (|¥T|) the corre- 
spondence i" — >• Pn{t) is a solution of eq. (gl]). 
Let us consider the dynamical system defined by differential equation 

(45) z'{t) — a2z{tY + aiz{t) + ao, 00,01,02 gM. 

It admits the general solution 

Vrtan{iVr(i + co)} 



(46) 



-Ol 



z(t) 



202 



where T = 4o2ao — a\, and Cg € M is fixed by the initial condition. By following 
the procedure described above, we obtain the discrete dynamical system 



(47) 

with solution 

(48) 



Zn+1 — Zn = 02 / ^ , ,, , 2fei^fc2 



/Cl,fc2=0 



fcilfca! 



{n — ki — k2)\ 



+ aiZn + Oo, 






fc=0 



{n-ky. 



where {/3k}keN is a number sequence, whose generating function is 
I fjk \/rtan| i\/r(a; + Co) ? - oi 



(49) 



Pk 



1 



fc! rfa;*^ 2o2 

Its regularized version of Borel type defines the map 

n 

(50) (n + l)m„+i-«;„ = 02 ^' (-l)'-^+'=^+"wfci«;fc, 

fcl,fc2=0 

whose solution is given by the finite Borel-type transform of z„: 



(n- fci - fc2)! 



-OlWn' 



Oo 



(51) 



Wn^Y. 



fe=0 



(n-fc)!' 



Comment 2. Besides the Borel regularization adopted in this work, other regu- 
larization procedures are possible (for instance, the Mittag-Lefher one). The other 
procedures a priori could lead to discrete dynamical systems different than (|44p or 

dsni). 
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4.3. A more general family of dynamical systems and associated recur- 
rences. Let us consider the dynamical system 

(52) —z = aNZ^, N eN gn eR. 

where z:M— >]RisaC°° function. The third system of the hierarchy is the 
dynamical system 

(53) Az = a^z*^. 

Its realization in the space (J^, +, •, *a+) can be performed in a similar way. There- 
fore, eq. (|53|) becomes on the lattice 



(54) Zn+i ~ Zn = a3 y^ 



f-^\ki+k2+k3+n 



~ ^k\ ^k2 '^k^ ' 



ki,k2,ks—0 

This system, as before, can be written in terms of the transformed variables (|26p . 
It has the form of a recurrence relation 

(55) (; + l)z;+i=a3 ^ zi^zi^zi_i^_i^. 

'l,!2=« 

;i+/2<i 
The solution of eq. (f52|) . for iV = 3, is provided (up to a sign) by 

(56) zit) = ^ 



\/2^/-a3t + Co ' 

where cq G K is an arbitrary constant fixed by the initial condition. Since our 
discretization preserves integrability, we get the following series solution of the 
dynamical system ([M|: 

^K'7\ sr^ Ik n\ a" 

(57) Zn = y —;=-, -—7 ,„, , ,x ,„ , n = 0,1,2,... 

where the first terms of the sequence {'^k}keTi are 

1 3 5 

70 = 1,71 == 2 '^2 = g,73 = Y^, 

35 63 

(58) j, = —, 75 = ^, etc. 

The generating function of this sequence is given by 

(59) -,, = 1^" 



dx'^ y/1 — X 

The procedure can be easily generalized: to each value of A^ G N, it corresponds 
naturally a specific sequence of rational numbers. 

The Borel regularization procedure, described above in the simpler case of qua- 
dratic dynamical systems, can be extended in a completely analogous way to the 
case (|54|) or to the maps arising from more general polynomial vector fields. The 
number sequence {6fc}fegN appearing in the explicit solution of the maps and of 
their Borel-transformed analogs are preserved. 

Consider now a general recurrence relation of the form 
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(60) {l + l)zi+i^aN ^ Z/1Z/2 •••Z;„Z/_ij_;2_...;_Ar. 

h+...+lK<l 

In order to find an exact solution of it, for every N , we interpret the recurrence 
()60p as tlie difference equation defining an abstract dynamical system, in the aux- 
iliary space A. Then, a particular solution of the recurrence can be obtained by 
associating with it the discrete map 



^ , (_l)ki-\ |-feiv+» 

Zn+l — Zn — ClN / : — j Tj ; ^fci 

•^ ' fci • • • fc^r 



(61) -Zfcj, 



fci!---!fc 



at: 



{n - ki - k2 - ■ ■ ■ kN)V 

obtained by means of the inverse transform (|26p . Then, we come back to the 
continuous system (|52[) . whose general solution has the form 

(62) z{t) = [{l-N)iaNt + co)]'^^'~''K 

Once we expand this solution and discretize it on the lattice C by following the 
procedure described above, we deduce a particular solution of the map (j6ip . which 
depends on the choice of the initial condition of ((5^ . in terms of a sequence of 
rational numbers; from (j26p we construct a solution of the recurrence (j60p . 
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